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We describe searches for B° meson decays to the charmless final states -k°K®K%, rjKgKg, and 
n'K° s K° s . The data sample corresponds to 467 x 10 6 BB pairs produced in e + e~ annihilation 
and collected with the BABAR detector at the SLAC National Accelerator Laboratory. We find no 
significant signals and determine the 90% confidence level upper limits on the branching fractions, 
in units of 10 -7 , B(B° -> ir°K%K%) < 9, B(B° -> r)K°Ks) < 10, and B(B° -» n' K° S K° S ) < 20. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 
The observation of mixing- induced CP violation in B° — > J/tjjKg decays as well as in the charmless 
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penguin-diagram dominated B° — > rj'K decays and 
of direct CP violation both in the neutral kaon system 
and in B° — > _ftT + 7r~ decays are in agreement with 
predictions of the standard model (SM) of electroweak 
interactions [Bj]. Further information about CP viola- 
tion and hadronic B decays can be provided by the mea- 
surement of branching fractions and time-dependent CP 
asymmetries in B decays to three-body final states con- 
taining two identical neutral spin zero particles and an- 
other CP eigenstate spin zero particle |6J. CP violating 
asymmetries have already been measured in B° decays to 
K°K°K° and to t^tt ^ 0, and a search has been 
performed in B — > rj'rj'K [9j. Other examples, in which 
study of time-dependent CP violation asymmetry might 
be particularly interesting, are the B° decays to tt K®K°, 
rjK^K®, and rj'KgK®. There are no theoretical estima- 
tions for the branching fractions of these SM-suppressed 
decay modes. Contributions from physics beyond the SM 
may appear in these decays. 

Among B meson decays to final states containing 
two kaons and an additional light meson, only B + — > 
K + K~it + has been observed, with a branching fraction 



of (5.0 ± 0.5 ± 0.5) x 10" B [10]. In this analysis an un- 
expected peak was observed around 1.5GeV/c 2 in the 
K + K~ invariant-mass spectrum. Studies of decays with 
two neutral or charged kaons in the final state, such as 
those presented herein, may help to elucidate the nature 
of this structure 



III. 



We present the results of searches for neutral B decays 
to charmless final states n°K°K°, rjK° s K° s and rj'K° s K% 
which arc studied for the first time. The results are based 
on data collected with the BABAR detector 121 at the 



PEP-II asymmetric-energy e + e~ collider located at the 
SLAC National Accelerator Laboratory. We use an inte- 
grated luminosity of 426 fb _1 , corresponding to 467 x 10 6 
BB pairs, recorded at the T(4S) resonance (center-of- 
mass energy y/s = 10.58 GeV) and, for the study of the 
background, 44 fb _1 collected 40 MeV below the reso- 
nance (off-peak). 

Charged particles from the e + e~ interactions are de- 
tected, and their momenta measured, by a combination 
of five layers of double-sided silicon microstrip detectors 
and a 40-layer drift chamber. Both systems operate in 
the 1.5 T magnetic field of a superconducting solenoid. 
Photons and electrons are identified with a CsI(Tl) crys- 
tal electromagnetic calorimeter. Charged particle identi- 
fication is provided by the average energy loss (dE/dx) in 
the tracking devices and by an internally reflecting, ring- 
imaging Cherenkov detector covering the central region 
(DIRC). A K/ir separation of better than four standard 
deviations (a) is achieved for momenta below 3 GeV/c. 
Detector details may be found elsewhere 12 1. 

The B daughter candidates are reconstructed through 
their dominant decays: rj — > 77 (^ 77 ), rj — ► 7r + 7r~7r (773^) 

WprJ where V -> 77, and 
h 7r~. We require the labo- 



where 7r° — > 77, rj 
rj' — > j o°7 (rj' pl ) where p° 



7/7T+7T 



ratory energy of the photons to be greater than 30 MeV 
for 7T° in 77 3Tr , 50 MeV for 7j 77 in rj' V7I ^, and 100 MeV 
for rf , and for 7r° and 7j 77 , produced directly from the 
B decay. We impose the following requirements on the 
invariant mass (in MeV/c 2 ) of the candidate final states: 
120 < m( 77 ) < 150 for tt°, 510 < 771(77) < 585 for ?i 77 
produced directly from the B decay, 490 < 771(77) < 600 
for 77 77 in rf^n, 538 < m(7r + TT~ n ) < 558 for rj 37T , 945 < 
m(7r+TT-rj) < 970 for 77^, 930 < m^ir-^) < 980 for 
77' , and 470 < to(7t + 7t~) < 980 for p°. Tracks from rj 
and rj' candidate decays are rejected if their particle iden- 
tification signatures from the DIRC and dE / dx are con- 
sistent with those of protons, kaons, or electrons. Can- 
didate Kg decays are formed from pairs of oppositely 
charged tracks with 486 < to(7t + 7t~) < 510 MeV/c 2 , 
a decay vertex x 2 probability larger than 0.001, and a 
reconstructed decay length greater than three times its 
uncertainty. 

We reconstruct the B meson candidate by combining 
two K® candidates and a ir°, rj, or 77' candidate. From the 
kinematics of the T(4S) decays we determine the energy- 
substituted mass toes = 



s — p B and the energy dif- 
ference AE = Eb — \\/~s, where {Eb, Pb) is the B meson 
4-momentum vector, and all values are expressed in the 
T(4S) rest frame. The resolution is 3.0 MeV/c 2 for toes 
and in the range (12-32) MeV for AE, depending on the 
decay mode. We require 5.25 < toes < 5.29 GeV/ c 2 and 
|A.E| < 0.2 GeV. 

Backgrounds arise primarily from continuum e + e~ — > 
qq events (q = u, d, s, c). We reduce these with a require- 
ment on the angle #t between the thrust axis of the B 
candidate in the Y(4S) rest frame and that of the rest 
of the charged tracks and neutral calorimeter clusters in 
the event [131 ]. The distribution is sharply peaked near 
I cos6*t| = 1 for qq jet pairs and is nearly uniform for B 
meson decays. The requirement is |cos#t| < 0.9. For 
the p° decays we also use | cos 6 P | where the helicity an- 
gle 9 P is defined as the angle between the momenta of a 
daughter pion and the 77', measured in the p° meson rest 
frame. For ?7 77 decays we use | cos 6 n \ where the decay 
angle 9^ is defined as the angle between the momenta 
of the most energetic daughter photon and the B° me- 
son, measured in the 77 meson rest frame. We require 
I cos6 p r v \\ < 0.9. Events are retained only if they con- 
tain at least one charged track in the decay products of 
the other B meson (B tag ) from the T(AS) decay. This 
requirement improves the precision of the determination 
of Btag thrust axis. The B° — * tt^K^K® decay has back- 
ground from B° -> TfiK®, with -> ir°K° s , which has 
the same final state as the signal mode. In order to sup- 
press this background, we define m(ir°K2.) as the closer of 
the two invariant mass combinations to the nominal D° 
mass [22j]. By requiring m(jT Kg) to be outside the range 
1.815-1.899 GeV/c 2 , we veto 80% of this background. 
We obtain the signal event yields from unbinned ex- 
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tended maximum likelihood (ML) fits. The observables 
used in the fit are AE, TOeSj and a Fisher discriminant 
T . The Fisher discriminant T [3| is a linear combina- 
tion of four event shape variables and |T|, the absolute 
value of the continuous output of a flavor tagging algo- 
rithm The event shape variables used for T are: 
the angles, with respect to the beam axis, of the B mo- 
mentum and the B thrust axis in the T(45*) frame, and 
the zeroth and second angular moments, Lq,2> of the en- 
ergy flow about the B thrust axis [III]. The moments 
are defined by Lj —J2iPi x |cos0i| J , where 9i is the an- 
gle, with respect to the B thrust axis, of track or neutral 
cluster i, and Pi is its momentum. The sum excludes the 
B candidate daughters. We use a neural network based 
technique to determine the flavor at decay of the 

Stag- 

The coefficients of T are chosen to maximize the sep- 
aration between the signal and the continuum back- 
ground. They are determined from studies of Monte 
Carlo (MC) [17| simulated signal data and off-peak data. 
Signal MC events are distributed uniformly across the 
Dalitz plot. Correlations among the ML input observ- 
ables are below 10%. The average number of candidates 
found per selected event is between 1.13 and 1.22, de- 
pending on the final state. We choose the candidate with 
the highest B vertex x 2 probability determined from a 
vertex fit that includes both charged and neutral parti- 
cles 18|. From simulated events we find that this algo- 
rithm selects the correct candidate in (92-98)% of the 
events containing multiple candidates, depending on the 
final state, and introduces negligible bias. 

We use a MC simulation to estimate backgrounds from 
other B decays, including final states with and without 
charm. These contributions are negligible for the r}'^^ 
mode. In all the other modes we introduce a non-peaking 
BB component in the fit. In the ir°K®K® analysis we also 
introduce a BB background component that peaks in 
toes and AE, to take into account the main contribution 
to background from B° — > K®KgK® decay mode. We 
consider three components in the likelihood fit: signal, 
continuum, and BB background. We have studied the 
possibility of misreconstruction of our B candidates. We 
divide signal events into two sub-components: correctly 
reconstructed (COR) signal and self cross-feed (SCF) sig- 
nal, where at least one B candidate daughter has been 
exchanged with a particle from the rest of the event. The 
signal component is split according to this classification. 
The fractions of SCF events are fixed in the fit to the 
values found in MC simulated events, which are in the 
range (10-21)%, depending on the final state. For the 
ifiK^K® decay mode, which has the lowest SCF fraction 
(6.6%), we use one signal component, comprising COR 
and SCF events. 

For each event i and component j, we define the prob- 



ability density function (PDF) 

V)=V ] {m^)V 3 {AV)V ] {F 
and the likelihood function: 



C 



JV 

n 

i=l 



(i) 



(2) 



where N is the number of reconstructed events and rij is 
the number of events in component j which is returned 
by the fit. We determine the PDF parameters from MC 
simulation of the signal and BB backgrounds, while we 
use toes and AE sideband data (5.25 < toes < 5.27 
GeV/c 2 , 0.1 < \AE\ < 0.2 GeV) to model the PDFs of 
continuum background. 

We parameterize "P(toes) as a Crystal Ball func- 
tion [l9( for the COR and SCF signal sub-components, 
an ARGUS function [20l ] for continuum and non-peaking 
BB background components, and by an ARGUS function 
plus an asymmetric Gaussian distribution for peaking 
BB background. The V(AE) distribution is described by 
an asymmetric Gaussian distribution plus an exponential 
tail (AGT) [2l[ for the COR signal sub-component, an 
asymmetric Gaussian distribution plus a linear Cheby- 
shev polynomial or an AGT for the SCF, and Chebyshev 
polynomials for continuum and BB background compo- 
nents. The distribution of T is described with an asym- 
metric Gaussian distribution plus a Gaussian distribution 
for the COR signal sub-component, an AGT function for 
SCF signal events, an asymmetric Gaussian distribution 
plus a linear Chebyshev polynomial for continuum, and 
an asymmetric Gaussian distribution for BB background 
sub-components. 

We allow the continuum-background PDF parame- 
ters to float in the fit. Large control samples of 
B~ — > £>°(if l?7r + 7r~7r )7r - decays are used to verify the 
simulated AE and toes resolution. Any bias in the 
fit, which mainly arises from neglecting the correlations 
among the discriminating variables used in the likelihood 
function definition, is determined from a large set of sim- 
ulated experiments. For each experiment, the qq back- 
ground and non-peaking BB background are drawn from 
the PDFs, and we embed the expected number of peak- 
ing BB background and signal events taken randomly 
from fully simulated MC samples. 

In TableUwe show, for each decay mode, the fitted sig- 
nal yields and their fit biases in numbers of events, the 
detection efficiencies, the product of daughter branching 
fractions, the significance S, and the measured branch- 
ing fractions. The detection efficiency is determined as 
the ratio of selected events in simulation to the num- 
ber generated. The significance is given in units of a. 
We determine the corrected signal yields from the fit- 
ted signal yields and their fit biases, estimated using 
simulations. We use these values, detection efficiencies, 
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TABLE I: Fitted signal yield in events and fit bias in events (ev), detection efficiency e (%), daughter branching fraction 
product significance 5 and measured branching fraction B with statistical error for each decay mode. For the combined 

measurements (in bold) we give S (with systematic uncertainties included) and the branching fraction with statistical and 
systematic uncertainties with the 90% CL upper limit in parentheses. 



Mode 


Yield (ev) 


Fit bias (ev) 


e(%) 


na (%) 




B(10- 7 


) 


n°K s K° s 


i -I 7 +16.2 
11 - 1 -14.5 


+ 1.0 ±0.7 


17.5 


47.9 


0.7 


2.7+37 ± °- 6 


« 9) 


rjK s K° s 


o r, +9.0 

-7.2 
n 9 +5.5 
-3.6 


+1.1 ±0.7 
+0.2 + 0.6 


17.5 
12.0 


18.8 
10.9 


0.3 
0.5 
0.5 


-4.7 
q q +9.0 
°-° -5.9 

2.1 it; ±1.2 


« io) 


rl K° K° 
W K°K° 

'lpj"-S-"-S 

T] K S K S 


2.4 ±tl 
13.4 tlti 


+0.1 + 0.4 
+4.7+1.1 


12.6 
15.9 


8.4 
14.1 


0.6 
0.6 
0.8 


4.6+1* 
8-3 tlH 
5.7t s e ° 5 ±3.4 


« 20) 



daughter branching fractions, and number of produced 
B mesons, assuming equal production rates of charged 
and neutral B meson pairs, to compute the branching 
fractions. The statistical error on the signal yield is the 
change in the central value when the quantity — 21n£ 
increases by one unit from its minimum value. The sig- 
nificance is the square root of the difference between 
the value of — 21n£ (with systematic uncertainties in- 
cluded) for zero corrected signal yield and the value at 
its minimum. We combine results from different sub- 
decay modes by adding the values of — 2 ln£. In order to 
account properly for systematic uncertainties when com- 
bining results from different sub-decays, we convolve the 
C of each sub-decay mode with a Gaussian distribution 
with mean equal to zero and width equal to the uncorre- 
cted systematic uncertainty of that decay mode. For the 
combined measurements we report the branching frac- 
tions, the statistical significances and the 90% confidence 
level (CL) upper limits. The 90% CL upper limit is taken 
to be the branching fraction below which lies 90% of the 
total likelihood integral in the positive branching fraction 
region. 

Figure [1] shows projections of 7r° K®K®, r)K®K®, and 
■q' K^K® candidates onto tties and AE for the subset 
of candidates for which the signal likelihood (computed 
without the variable plotted) exceeds a mode-dependent 
threshold. 

The main sources of systematic error include uncer- 
tainties in the detection efficiencies, the PDF parameters, 
and the maximum likelihood fit bias. We assign system- 
atic uncertainties (13-20%) on the detection efficiencies 
due to non-uniformity of the efficiencies over the Dalitz 
plot. This contribution is taken to be the ratio between 
the standard deviation of the efficiency distribution over 
the Dalitz Plot to its mean value. For the signal, the 
uncertainties in the PDF parameters are estimated by 
comparing MC and data control samples. Varying the 
signal PDF parameters within these uncertainties, we es- 
timate the yield uncertainties of 0-2 events, depending 
on the mode. The uncertainty from the fit bias is taken 
as the sum in quadrature of one-half the correction (1-3 
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FIG. 1: B° candidate tties and AE projections for 
Tr°K^Kl(a,b), nK° s K° s (c,d), and for n' K° S K° S (e,f) with the 
sub-decay modes combined. Points with errors represent the 
data, solid curves the full fit functions and dashed curves the 
background functions. These plots are made with a require- 
ment on the likelihood in order to enhance signal to back- 
ground ratio. 



events) plus the statistical uncertainty on the correction 
itself. We assign a systematic error of 0.1-0.4 events, de- 
pending on the mode, due to non-uniformity of the SCF 
fraction over the Dalitz plot. Uncertainties of the effi- 
ciency found from auxiliary studies include 0.8% x N t 
where Nt is the number of tracks in the B candidate. A 
systematic uncertainty of 1.8% and 3.0% is assigned to 
the single photon and 7r°/77 77 meson reconstruction effi- 
ciencies, respectively. There is a systematic error of 0.9% 
for the reconstruction efficiency of each K®. The uncer- 
tainty on the total number of BB pairs in the data sam- 
ple is 1.1%. Uncertainties on the B daughter branching- 
fraction products (3.5-4.9)% are taken from Ref. [22}. 
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In conclusion we have searched for the B decay modes 
to Tr°KjK°, T]K°K° S and T]'K S K S with a sample of 467 x 
10 6 BB pairs. We find no significant signals and set 90% 
CL upper limits for the branching fractions: B(B° — > 
tt°K°K° ) < 9 x 1(T 7 , B(B° -> T)K°K°) < 10 x 10~ 7 , 
and B(B° -> j/if°iT° ) < 20 x lO" 7 . 
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